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Abstract 

Parity-odd domains, corresponding to non-trivial topological solutions of the QCD vacuum, 
might be created in relativistic heavy ions collisions. These domains are predicted to lead to 
charge separation along the system orbital momentum of the system created in non-central col- 
lisions. Three-particle mixed harmonics azimuthal correlator is a !P even observable but directly 
sensitive to the charge separation effect. Using this observable to analyze Au-hAu and Cu+Cu 
collisions at yfs]^ = 200 and 62 GeV, STAR detects a signal consistent with several of the theo- 
retical expectations. Possible contributions from effects not related to parity violation are studied 
with existing event generators, which fail to describe the data. Future directions in studying the 
effect are discussed. 



1. Introduction. Quantum Chromodynamics (QCD) is the theory of strong interactions. Per- 
turbative QCD is firmly established and thoroughly tested experimentally. In the non-perturbative 
sector, QCD links chiral symmetry breaking and the origin of hadron masses to the existence of 
topologically non-trivial classical solutions describing the transitions between the vacuum states 
with different Chern-Simons numbers. Quark interactions with topologically non-trivial clas- 
sical gluonic fields change the quark helicity and are !P and CP odd. It was suggested in HI 
that metastable !P and CP odd domains, characterized by non-zero topological charge, might be 
created in ultra-relativistic heavy ion colUsions. The possibility for an experimental detection of 
this local strong parity violation was discussed in ||Tl|2l[3l. More recently, it was noticed ||4l|5l 
that in non-central collisions such domains can demonstrate themselves via the asymmetry in the 
emission of charged particle w.r.t. the system's angular momentum. Such charge separation is 
a consequence of the difference in the number of particles with positive and negative helicities 
positioned in the strong magnetic field (~ 10'^ T) of a non-central nuclear collision It4ji6j, the 
so-called chiral magnetic effect. The same phenomenon can also be described in terms of the 
induced electric field that is parallel to the static external magnetic field, which occurs in the 
presence of topologically non-trivial vacuum solutions fT). 

Since the direction of the separation may vary event by event in accord with the changing 
sign of the topological charge of the domain, the observation of the effect is possible only by 
correlation techniques. Such an observable, !P-even, but directly sensitive to the charge sepa- 
ration effect, has been proposed in |8| and is based on 3-particle mixed harmonics azimuthal 
correlations. 

Phenomenologically, the charge separation can be described by adding a !P-odd sine term to 
the Fourier decomposition of the particle azimuthal distribution with respect to the reaction plane 
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angle, W^p, which is often used in the description of anisotropic flow |[8l : 



— - — oc 1 + 2vi cos(A(f>) + 2v2 cos(2z)0) + ... + 2a± sin(zJ0) + ... , 



(1) 



where Acp - (4> - ^rp) is the particle azimuth relative to the reaction plane, vi and V2 account 
for directed and elliptic flow. Parameters a_ = -a+ describe the !P-violating effect. The sign 
of fl-t varies event to event following the fluctuations in the domain's topological charge, and 
on average is zero, (a-t) = 0. The observation of the effect is possible via correlations, e.g. 
(ciaa^, where a and [5 denote the particle type. To measure ^a^a^^, it was proposed [SJ to use 
the correlator: 



This correlator represents the difference between correlations projected onto an axis in the reac- 
tion plane and the correlations projected onto an axis perpendicular to the reaction plane. The 
key advantage of using such a difference is that it removes all the correlations among particles 
a and y6 that are not related to the reaction plane orientation |i9, J_OJ . Only the parts of such 
correlations that depend on azimuthal orientation with respect to the reaction plane remain as 
backgrounds, denoted as [Bom - Bi„]. Note that the contribution given by the term ^vi „vi^^^ can 
be neglected because directed flow averages to zero in a rapidity region symmetric with respect 
to mid-rapidity, as used in this analysis. 



According to Refs. |l4l|5]|6| one expects the following features of the correlator 



• Magnitude: Estimates ID predicted a signal \a\ ~ Q/N^*, where Q = 0, ±l,+2, ... is the 
topological charge and A^;^* is the positive pion multiplicity in one unit of rapidity - the 
typical scale of such correlations. More accurate estimates |6| including the strength of 
the magnetic field and topological domains production rates, were found to be close to the 
same number, of the order of 10"^ for mid-central collisions at top RHIC energies, which 
corresponds to 10""* for the correlator (aaapj. 



• Charge combinations: In the absence of medium effects, one expects {a^a+) - {a^aJ) - 
— {a+aJ) > 0. In a dense medium suppression of the back-to-back correlations may oc- 
cur [6|: {a+a+) = (a^a^) » - {a^aJ). The effect of strong radial flow can further modify 
this relation such that the opposite charge correlations can become even positive. 

• Centrality dependence: The correlator is expected to follow a 1/N dependence (typical 
for any kind of correlations due to clusters; is the multiplicity) multiplied by a factor 
accounting for the variation of the magnetic field. Thus at large centralities the effect 
should decrease with centrality somewhat faster than 1 /N. 

• Rapidity dependence: The correlated particles come from a domain of the order of 1 fm, 
and it is expected that the correlations would have a typical hadronic width in At] - \Tia-t]js\ 
of the order unity. 

• Transverse momentum dependence: Local parity violation is non-perturbative in nature 
and one would expect the main contribution to the signal at pt < 1 GeV/c, but the actual 
limits might be affected by the radial flow. 
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• Beam species dependence: The effect should be proportional to the square of the nu- 
clear charge, but the atomic number A dependence is not well understood. The suppression 
of the back-to-back correlations should be smaller in collisions of lighter nuclei. 

• Collision energy dependence: The effect might be stronger at lower energies, as the time 
integral of the magnetic field is larger At the same time, the charge separation effect is 
expected to depend strongly on deconfinement and chiral symmetry restoration |6|, and 
the signal might be greatly suppressed or completely absent at an energy below that at 
which a quark-gluon plasma can be formed. 

2. Data. The results are based on 14.7M Au-hAu and 13.9M CuH-Cu events at y[s^-2QQ GeV, 
and 2.4M Au-hAu and 6.3M Cu+Cu events at ^Js^-62 GeV recorded with the STAR detec- 
tor ifTTI at RHIC. Charged particle tracks were reconstructed in a Time Projection Chamber 
(TPC), \ri\ < 1.0, operated in a solenoidal magnetic field of 0.5 T. The main TPC is supplemented 
by two Forward TPCs, which cover pseudorapidity intervals 2.7 < \ri\ < 3.9. A minimum bias 
trigger was used during data-taking. The centrality of the collision is determined according to the 
recorded multiplicity of charged particles in \t]\ < 0.5. The correlations are reported for charged 
particle in the region I;;! < 1.0 with p, > 0.15 GeV/c. The results integrated over transverse 
momentum have a cut of p, <2 GeV/c. 

3. Method. In the three-particle correlation technique, used in this analysis, the role of the 
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Figure 1 : A comparison of the correlations obtained 
by selecting the third particle from the main TPC or 
from the Forward TPCs after the results scaled by the 
flow of the third particle. The shaded areas represent 
the uncertainty from V2,c scaling (see text for details). 
The error bars are statistical. 



event plane is played by the third particle that enters the correlator with double the azimuth IfTOl 
\T2\ . Under the assumption that particle c is correlated with particles a and /3 only via common 
correlation to the reaction plane, one has: 

<COS(0„ +<pi3- 24>c)) = <COS(0fl +(f>i3- 2Wrp)) V2,c, (4) 

where V2,c is the elliptic flow value of the particle c. We check this assumption by using particles 
c from different detectors that exhibit different elliptic flow. Figure[T]compares the three-particle 
correlator divided by V2,e, when the third particle is selected from the main TPC with those when 
it is selected from the Forward TPCs, where elliptic flow is significantly smaller 1 13 1. The shaded 
bands in Fig. [T] and below illustrate the systematic error due to uncertainty in V2 measurements. 
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The upper (in magnitude) limit is obtained with V2{4) and the lower limit from V2{2), the mid- 
point is calculated using V2{FTPC} Lii J4|- Whenever V2{4) values are not available the upper 
limits are obtained assuming that the measurements with FTPC suppress only 50% of the non- 
flow contribution. 

A very good agreement between the same charge correlations in Fig.[T]justifies the assump- 
tion Eq. |4] But the opposite charge correlations are too small in magnitude to conclude on 
validity of the factorization. Similarly, in the most peripheral collisions, the statistical errors are 
large, which also prohibits making a definite conclusion. In these cases, the factorization can be 
broken due to contribution of three-particle clusters, which we estimate with the help of event 
generators. We proceed below assuming (cos{<pa + (p/3 - 2Wrp)^ - (cos((pa + <P/3 ~ '^'Pc)) h'l.c but 
indicate in all plots the HIJING [ISJ direct three -particle correlation results. The latter can be 
considered as an estimate of the systematic uncertainty from correlations not related to the re- 
action plane. (Using the event generator UrQMD yields both opposite-charge and same-charge 
correlations at least a factor of two lower than those predicted by HIJING.) Note that in principle 
such uncertainty can be suppressed by taking particle c farther in rapidity from the particles a 
and/?. 

4. Results. Figure ji^ presents the correlator (cos(if>a + (pp- IWrp)^ for Au-hAu and CuH-Cu 
collisions at yiJviv=200 GeV. {+, +) and (-, -) correlations are found to be the same within 
statistical errors and are combined together as same-charge correlations. Opposite-charge corre- 
lations are relatively smaller than same-charge correlations. This observation led to the proposal 
that back-to-back correlations may be suppressed due to the opacity of the medium as discussed 
in the introduction. The difference in magnitude between same and opposite sign correlations 
is considerably smaller in Cu+Cu than in Au+Au, qualitatively in agreement with the scenario 
of stronger suppression of the back-to-back correlations in Au-i-Au collisions. In Fig. [2^ and 
below, error bars show the statistical errors. Note again that presenting the results in this section 
we assume the factorization of correlator, Eq. |4] The possible error due to this assumption is 
denoted by the thick lines in Fig.|2]and subsequent figures. Figure|2|5 shows results for collisions 
at yjJviv=62.4 GeV. The signal is similar in magnitude, with slightly more pronounced opposite- 
charge correlations compared to those at y/J^N-lOO GeV. This is consistent with weaker sup- 
pression of opposite-charge correlations in the less dense 62 GeV system. 

The correlations are weaker in more central collisions compared to more peripheral colli- 
sions, which partially can be attributed to dilution of correlations which occurs in the case of 
particle production from multiple sources. The somewhat stronger correlations in CuH-Cu colli- 
sions than in Au+Au for the same centrality also may have similar explanation. To compensate 
for the dilution effect we show in Fig. |3] results multiplied by the number of participants. The 
decrease of the correlations in most central collisions is expected as the magnetic field weakens. 
The same and opposite sign correlations clearly exhibit very different behavior. The opposite 
sign correlations in Au+Au and Cu+Cu collisions are found to be very close at similar values of 
Npart in rough qualitative agreement with the picture in which their values are mostly determined 
by the suppression of back-to-back correlations. 

Figure shows the dependence of the signal on the diff'erence in pseudorapidities of two 
particles, At] - \ria - rj^l, for 30-50% centrality. The signal has a typical hadronic width of about 
one unit of pseudorapidity. Figure |^ shows the dependence of the signal on the sum of the 
transverse momentum (magnitudes) of the two particles for these same centralities. The signal 
is not concentrated in the low p, region as naively might be expected for !P-violation effects. We 
find also that the correlation depends very weakly on \p,^a - Ptfi\ (not shown). This excludes 
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Figure 2: ^cos(0a + (pp - 2¥rp)^ in Au+Au and Cu+Cu collisions at (a) ^iivyv = 200 GeV and (b) -\jSNN = 62 GeV. 
Thick solid (Au+Au) and dashed (Cu+Cu) lines represent possible non-reaction-plane dependent contribution froin 
many-particle clusters as estimated by HIJING. 
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Figure 3: Au+Au and Cu+Cu collisions at 
yiAiAi=200 GeV. The correlations are scaled with the 
number of participants and shown as a function num- 
ber of participants. Thick lines as in Fig.|2] 
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Figure 4: Au+Au at 200 GeV. The correlations dependence (a) on pseudorapidity separation At] = \t]„ — r]p\, and (b) on 
sum of the magnitude of transverse momenta. 



quantum interference or Coulomb effects as possible explanations for the signal. There are no 
specific theoretical predictions on this dependence, though naively one would expect that the 
signal should not extend to large values of - p,fi\. 

5. Physics background. The correlator ^cos(0„ + 0/3-2 ^'sp)^ is a !P-even observable and 
can exhibit a non-zero signal for effects not related to !P-violation. Among those are processes 
in which particles a and /3 are products of a cluster (e.g. resonance, jet, di-jets) decay, and the 
cluster itself exhibits elliptic flow |10| or decays (fragments) differently when emitted in-plane 
compared to out-of-plane. If "flowing clusters" are the only contribution to the correlator, one 
can write: 

NchnLNp^ 

dust ^2,c.lust-> (5) 

i V poil s 
evem 

where {■■■)ciust indicates that the average is performed only over pairs consisting of two daughters 
from the same cluster Estimates of the contribution of "flowing resonances", based on Eq.|5]and 
reasonable values of resonance abundances and values of elliptic flow, indicate that they should 
not produce a fake signal. 

To study the background contribution in greater detail we used event generators MEVSIM |l76l. 
UrQMD ini and HIJING 1 15|. The results are presented in Fig.|5] MEVSIM includes only cor- 
relations due to elliptic flow and resonance (0, A, p, w, and K*) decay. HIJING is also run with an 
added "afterburner" to add elliptic flow as experimentally observed. No generator gives qualita- 
tive agreement with the data; the model values of ^cos(0q, + (pjj-l 'Frp)^ are significantly smaller 
in magnitude than what is seen in the data, and the correlations calculated in these models tend to 
be very similar for same and opposite sign correlations. In Fig. |5] UrQMD points are connected 
by dashed lines to illustrate that the "reference line" for strong parity correlations might be not 
at zero. 

Here, we also mention that such effects as directed flow fluctuations and global polariza- 
tion Il20ll2n should not have any significant contribution to the measurements. 
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Figure 5: Data comparison to simulation results for 
200 GeV Au+Au. Blue symbols mark opposite- 
charge correlations, and red are same-charge. Mark- 
ers connected by solid lines represent the data. Ac- 
ceptance cuts of 0.15 < Pi < 2 GeV/c and \t]\ < 1.0 
were used in all cases. In simulations the true reac- 
tion plane from the generated event was used. Thick 
solid lighter colored lines represent non reaction-plane 
dependent contribution as estimated by HIJING. Cor- 
responding estimates from UrQMD are about factor of 
two smaller. 



6. Future directions. Taking into account the importance of the question, one can envision 
a dedicated program for estabUshing the nature of the signal and further detail study. From the 
theoretical point of view, the calculation of the dependence on centrality and system system size 
looks fully doable though requires significant computing and man power (e.g. 3d hydrodynamics 
is needed for the calculation of the magnetic field). Detailed predictions on the transverse mo- 
mentum and particle type dependence of the effect also will be essential in differentiating it from 
possible "background" contributions. A good theoretical understanding of the "background" 
correlations itself is also required, as at present all event generators lack a good description of 
correlation results. 

A number of future experiments and analyses are naturally suggested by STAR results. One 
of them is the dependence of the signal on the energy of the colliding ions, which can be ad- 
dressed, for example, during the RHIC beam energy scan The charge separation effect is ex- 
pected to depend strongly on the formation of a quark-gluon plasma |6|, and the signal might be 
greatly suppressed or completely absent at an energy below that at which a quark-gluon plasma 
can be formed. Identified and multiparticle correlations studies also will be available with larger 
statistics. Those will be important to test a specific predictions, e.g. such as topological cluster 
decays in equal number of <7^-pairs of all flavors. The charge separation dependence on the mag- 
netic field Q can be tested with collision of isobaric nuclei, such as I^Ru and ^^Zr that were used 
at GSI (221. Collision of isobaric nuclei will be also very interesting in relation to the directed 
flow studies (in the latter case one needs asymmetric collisions). 

Correlation measurements from RHIC [! 18 1 and earlier measurements at ISR (see review |fT9l ) 
indicate that cluster formation plays an important role in multiparticle production at high ener- 
gies. These clusters, with a size inferred in |18 | to be 2.5-3 charged particles per cluster, may 
account for production of a significant fraction of all particles. It will be interesting and impor- 
tant to establish how these clusters are related to the topological clusters as suggested in 1231 
and ||24 I ("turning points"). Note recent progress in describing of the soft Pomeron as a multi- 
instanton ladder ESl . which also suggests an important role played in multiparticle production 
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by topologically nontrivial gluonic configurations. 

7. Summary. Charge separation due to quark interaction with topologically non-trivial glu- 
onic configurations in a strong magnetic field of non-central heavy ion collision may provide a 
unique opportunity for a direct observation of the topological structure of QCD. The predictions 
are well within reach of the experiment. STAR has performed analysis of Au-i-Au and CuH-Cu 
collisions at -^5jviv=200 and 62 GeV using three-particle correlations that are directly sensitive to 
the local !P-violation effects in heavy-ion collisions. The results are reported for different particle 
charge combinations as a function of collision centrality, particle separation in pseudorapidity, 
and particle transverse momentum. Qualitatively the results agree with the magnitude and gross 
features of the theoretical predictions for local !P-violation in heavy-ion collisions, but the sig- 
nal persists to higher transverse momenta than expected [6 J. The particular observable used in 
the analysis is 7^-even and might be sensitive to non-parity-violating effects. None of the stud- 
ies we have performed so far have revealed a background source that can explain the observed 
same-sign correlations. 

Better theoretical calculations of the expected signal and potential physics backgrounds in 
high energy heavy-ion collisions are essential for confirmation and experimental study of this 
phenomenon. 

A very exciting future program dedicated to detail study of the efifect is emerging. 
Acknowledgments. Discussions with D. Kharzeev and E. Shuryak are gratefully acknowl- 
edged. 
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